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ABSTRACT 

We introduce the MUSCEL Program (MUltiwavelength observations of the Structure, 
Chemistry and Evolution of LSB galaxies), a project aimed at determining the star- 
formation histories of low surface brightness galaxies. MUSCEL utilizes ground-based 
optical spectra and space-based UV and IR photometry to fully constrain the star- 
formation histories of our targets with the aim of shedding light on the processes 
that led low surface brightness galaxies down a different evolutionary path from that 
followed by high surface brightness galaxies, such as our Milky Way. Here we present 
the spatially-resolved optical spectra of UGC 628, observed with the VIRUS-P IFU 
at the 2.7-m Harlen J. Smith Telescope at the McDonald Observatory, and utilize 
emission-line diagnostics to determine the rate and distribution of star formation as 
well as the gas-phase metallicity and metallicity gradient. We find highly clustered 
star formation throughout UGC 628, excluding the core regions, and a log(0/H) 
metallicity around —4.2, with more metal rich regions near the edges of the galactic 
disk. Based on the emission-line diagnostics alone, the current mode of star formation, 
slow and concentrated in the outer disk, appears to have dominated for quite some 
time, although there are clear signs of a much older stellar population formed in a 
more standard inside-out fashion. 

Key words: galaxies: evolution, galaxies: abundances, galaxies: spiral, dust, extinc¬ 
tion 


1 INTRODUCTION 


Low surface brightness (LSB) spiral galaxies represent a sig¬ 
nificant gap in our understanding of star formation: They 
are gas-rich disks which span the same range of baryonic 
masses as their high surface brightness (HSB) counterparts, 
yet, for reasons not fully understood, they have largely re¬ 
tained their atomic gas instead of producing stars. The dif¬ 
ficulty in detecting LSB galaxies makes their cosmological 
significance challenging to estimate; studies based on optical 
surveys ( McGaugh et al.|1995 O’Neil fe Bothun|2000 Tra- 


chternach et al. 2006 ) indicate that LSB galaxies account 


for up to 50% of the galaxy-bound baryons in the universe, 


fraction closer to 10%. 


while Hayward et al. (2005 I use supernovae to estimate that 


Attempts to place LSB galaxies within the context of 
hierarchical formation using broad-band colors have been 
hampered by their paradoxical properties, such as blue col¬ 
ors but faint stellar populations and low star-formation rates 
despite gas rich disks. The only universal conclusion that 


has been reached is that the star-formation histories of LSB 
galaxies differ from their HSB counterparts. 

Among the many descriptions of LSBs, a favored phe¬ 
nomenological explanation for the blue colors and low star- 
formation rates is that LSB galaxies follow a burst-and- 
quench cycle, where small regions sporadically erupt with 
short bursts of star formation just frequently enough to keep 
the galactic disk blue without building up a large stellar 


population (Zackrisson et al. 2005 Vorobyov et al. 2009). 


Observations of patchiness in Ha images of LSB galaxies 
(e.g., Auld et al.||2006 Kim|2007 1, consistent with the brief 
but concentrated star formation proposed by Vorobyo v et al.| 
12009), support this hypothesis. 


The burst-and-quench scenario is in line with a signif¬ 
icant clue from the ubiquitously low gas surface density in 
LSB galaxies (despite a globally high M ga , a /M* ratio), with 
p g as rarely climbing above the Kennicutt-Schmidt threshold 
(Kennicutt|1998|). Several authors (e.g., de Blok et al. 1996 


|Mihos et al |1997| |Schombert et al.|2013 1 argue that the gas 

density is too low to allow for self-gravitating global pat- 
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terns, leaving star formation (such as it is) to be governed by 
local fluctuations. For example, while 50% of HSB galaxies 


host nuclear star formation, Schombert et al. (2013) found 


that nuclear star formation is extremely rare in their sample 


of 54 LSB galaxies. Likewise, Mihos et al. (1997) point out 


that barred LSB spirals appear to make up 3-4% of several 
samples of LSB galaxies (McGaugh fc Bothun|1994 de Blok 


et al.|[l995 Impey et al.||1996 l, while 30% of isolated HSB 

spirals have bars (Elmcgrccn et al. 1990). 

The role of global patterns, however, is far from clear, as 


Kim (2007) found that in a sample of 13 spiral and/or barred 
LSB galaxies seven show centralized star formation and 
six show correlation between star formation and the arms 
and/or bars, suggesting that global patterns may be signifi¬ 
cant. Global patterns may also be more of a side effect than 
the ultimate cause; several authors (e.g., Bothun et al.|1993 


Mo et ak||T994) have hypothesized that LSB galaxies have 


low star-formation rates because they are more isolated, and 
interactions are a major source of galaxy-wide patterns. Sup¬ 


porting this idea, Rosenbaum & Bomans (2004) found that 


LSB galaxies in SDSS data are much more likely to be found 
on the edges of large scale structure, likely having formed in 
the voids. 

Another significant clue to the star-formation histories 
of LSB galaxies may be found in their dust and metal con¬ 
tent. With modest stellar populations and high gas fractions 
LSB galaxies are generally assumed to be largely unevolved 
and relatively free of dust and metals. This assumption is 
circumstantially supported by weak or absent thermal dust 


emission in the mid- and far-IR Hinz et al. (2007), however, 


since thermal dust is more steeply correlated with UV heat¬ 
ing than it is with actual dust content, the absence of ther¬ 
mal dust emission is likely linked to the low star-formation 
rate. Indeed, Hinz et al. ( |200 7) note that in their sample 
of five LSB galaxies thermal dust emission is found only in 
those with HII regions. 

In a comprehensive study [McGaugh (19941 spectro¬ 
scopically surveyed HII regions in 22 LSB galaxies, and 
found that the reddening and log(0/H) metallicities span 
a wide range, albeit with averages lower than HSB galaxies. 


Vorobyov et al. (20091 point out that the burst-and-quench 


scenario leaves the signature of metallicity variations across 
the galactic disk if LSB galaxies are truly unevolved; since 
no such fluctuations are seen in observations, LSB galaxies 
appear to have been forming stars for quite some time, albeit 
very slowly. 

To thoroughly address the star-formation histories of 
LSB galaxies, we have undertaken a multiwavelength study 
of a sample of gas-rich Milky Way mass LSB spirals cen¬ 
tered around optical IFU observations; this is the first in 
a series of papers which will closely examine current and 
past star formation in LSB galaxies with the aim of placing 
this difficult to observe galaxy class in the context of galaxy 
formation and shedding light on the physics which regulate 
star formation. 

In this paper we examine the optical emission-line 
strengths and ratios of HII regions in our first target, 
UGC 628. The information available via the emission lines 
allows us to place this LSB galaxy in the context of the 
various models and scenarios proposed, and will supplement 
the results in an upcoming series of papers. Specifically, we 
examine the optical emission-line properties of UGC 628, fo¬ 


cusing on the [OIIJA3727, H/3, [OIIIJA4959, and [OIII]A5007 
fines. These fines, particularly the H/3 line, allow us to ad¬ 
dress the rate and distribution of current star formation in 
UGC 628, while the fine ratios act as windows into the star- 
formation history in UGC 628. 

In Section [2] we describe our target selection, observa¬ 
tions, and basic data reduction, in Section[3]we describe the 
methodology used to process the reduced data into narrow- 
band images, in Section [4] we compare our findings to those 
in other studies of LSB and HSB galaxies, and in Section [5] 
we summarize our findings and place UGC 628 in the con¬ 
text of galaxy evolution. 


2 OBSERVATIONS 

Our program, MUSCEL (MUltiwavelength observations of 
the Structure, Chemistry and Evolution of LSB galaxies), 
is designed to unambiguously determine the star-formation 
histories of LSB galaxies, a goal which faces the signifi¬ 
cant challenge that the parameters of model galaxies are 
often degenerate over observations in any one wavelength 
regime. To combat this issue, we have designed our program 
in three segments: space-based UV observations (mostly via 
the Swift UVOT instrument), which fully constrain the cur¬ 
rent star-formation rate; archival Spitzer IRAC observa¬ 
tions, which fully constrain the integrated star-formation 
history; and ground-based optical spectra, which are sen¬ 
sitive to intermediate age indicators, such as the shape of 
the stellar continuum, Balmer absorption fines, and the 
4000A break. 

In order to match the ~2" angular resolution of our 
broad band UV and IR images, we chose to make our spec¬ 
tral observations with the VIRUS-P integral field unit (IFU), 
a multi-fiber spectrograph on the 2.7-m Harlen J. Smith 
Telescope at Mt. Locke, which has the capability to provide 
a spectrum at every location in the field of view, essentially 
giving us an image of our target at every wavelength within 
the range of our instrument. The VIRUS-P wavelength range 
encompasses several bright optical emission lines which, by 
themselves, provide significant clues to the star-formation 
history of our first target, UGC 628. An upcoming paper 
will examine the star-formation history fitted against the 
full UV/optical/IR SED, but here we present an analysis of 
the emission lines alone. 


2.1 Target Selection and Background 


We chose targets drawn from the LSB galaxies cataloged in 


Kim (2007) and McGaugh & Bothun (1994) that fit com¬ 
fortably within the 1.7 x 1.7 VIRUS-P field-of-view while still 
spanning several 2" seeing-limited resolution elements (typ¬ 
ical of Mt. Locke). To maximize galaxy coverage, minimize 
dust obscuration, and minimize the number of distinct stel¬ 
lar populations along each line-of-sight, we avoided galaxies 
that are edge-on (i > 85°). Additionally, since the primary 
objective of our program is an evaluation of candidate star- 
formation histories, we restricted our targets to LSB galax¬ 
ies with archival Spitzer IRAC observations to constrain the 
mature stellar population. The first target in our program 
with completed ground-based spectra is UGC 628, which is 
the focus of this paper. 
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10" 5 kpc 

Figure 1. (a, left) SDSS image of UGC 628. (b, right) An H/3 
image of UGC 628 overlayed with stellar continuum isophotes. 


With a central B-band surface brightness of /ro = 
23.1 mag arcsec -2 (Kim 2007) UGC 628 falls well within 
the LSB category, defined as fi o > 22 by McGaug h et al.| 
(1995]), though it is by no means an extreme member of 
the population (/.to > 24). Using broad-band photometry 


Kim (2007) estimate log(M*/M©) = 10.65 or 10.80, depend¬ 


ing on the IMF assumed, making it a near match to the 
log(M*/M©) = 10.7 for the Milky Way (|Flynn et al.||2006 
McMillan|20TT |. 


UGC 628 is an extreme late-type spiral (see Fig- 


couleurs et al. 

1991 

1 . With a redshift of 2 = 0.018166 

(Schneider et al. 

1990 

), the proper distance to UGC 628 is 


approximately 78Mpc. Like many LSB galaxies, UGC 628 
is dark matter dominated, with nearly solid-body rotation 
out to ~3 kpc ( |de Blok fe Bo sma 20 02)), but, unlik e most 
LSB disks, UGC 628 is a barred spiral. Mihos et al. (19971 
point out that only a few percent of LSB disks show bars, 
compared to the 30% of HSB spirals reported in the RC2 
catalog (El megreen et al.|19 90). 

The bar in UGC 628 is itself also extremely unusual in 
that it is a slow bar, meaning that bar pattern speed is slow 
enough that the corotation radius is more than 1.4 times the 
semi-bar length (Chemin & Her nandez| 2009). This is one of 
only three slow bars known, the others being a stellar bar in 
the dwarf irregular galaxy NGC 3741 (Banerjee et al. 2013) 
and a purely gaseous slow bar in the dwarf irregular galaxy 
NGC 2915 ( |Bureau et al.|1999| |. 


fold. First, because the VIRUS-P IFU has a 1 /$ fill factor, a 
minimum of three pointings is needed to completely fill the 
field-of-view. Second, since 4"16 VIRUS-P fibers undersam¬ 
ple the 2" seeing (typical for Mt. Locke), over dithering by 
a factor of two is necessary to approach the seeing-imposed 
limit (see Section 2.51. 

In addition to the VIRUS-P IFU frames, we took con¬ 
temporaneous frames with the fixed-offset guide camera. 
The guide camera data were used to calibrate the relative 
astrometry of the IFU data, which enabled reconstruction 
of the field-of-view (see Section 2.5). 


2.3 Basic Reduction 

In addition to the science frames we also collected bias, flat- 
field, and wavelength calibration data at the beginning and 
end of each night. We median-combined all of the bias frames 
taken on each of the four nights to create a nightly master 
bias frame, with which we removed the bias from each of 
the non-bias frames on each respective night. The remain¬ 
ing flat-fields, target frames, sky frames, and wavelength cal¬ 
ibration frames were then cleaned to eliminate cosmic rays. 

The cleaned wavelength comparison frames were 
median-combined to produce high S/N cosmic-ray-free mas¬ 
ter wavelength comparison frames for each evening and 
morning, which we then compared to a Hg-Cd line list to 
create an evening and morning wavelength solution for each 
night. 

Using this wavelength solution the cleaned flat, target, 
and sky frames were then wavelength-rectified and collapsed 
— that is, within each x column all of the pixels associated 
with each fiber were summed to a single value, and those 
values were binned with regular steps in wavelength to align 
the fibers. To account for variations in the telescope and 
spectrograph over the course of the night, we used the wave¬ 
length solution (evening vs. morning) taken closest in time 
to each respective frame. 

Similarly, to create a master flat-field for the evening 
and morning of each night, we normalized and median com¬ 
bined the collapsed twilight flat-fields. All 246 spectra in 
each sky flat were divided by the underlying sky spectrum, 
which we determined by taking the median value of the 246 
spectra. Each of the target and sky frames were flattened 
and collapsed using the flat field, evening or morning, that 
was closest in time. 


2.2 Data Acquisition 

We collected our data during a four night observing run in 
November 2013. The weather was mostly clear for one and 
one half of these nights, and during this time we completed 
observations of UGC 628. 

The VIRUS-P IFU consists of 246 fibers permanently 
packed into a hexagonal bundle and placed at the focal plane 
with the spectrally dispersed images of the fibers refocused 
onto a 2048x2048 CCD. Since the aim of our project is to 
determine the shape of the stellar continuum in our target 
galaxies, we used the VP1 grating with a 5.1 A resolution to 
maximize wavelength coverage, spanning roughly 3525A to 
5775A . 

We planned our observations in a six-point dither pat¬ 
tern. The purpose of the six-point dither pattern is two- 


2.4 Sky Subtraction 

The surface brightnesses of our program’s targets are all 
fainter than the background of even the darkest skies, mak¬ 
ing precise sky subtraction paramount. Fortunately, our tar¬ 
gets are all smaller than the VIRUS-P field-of-view, leaving 
the fibers at the edge of the field-of-view as contemporane¬ 
ous sky measurements. In our final data products, we adopt 
the median of the faintest forty fibers in each target ex¬ 
posure as the sky background spectrum, and subtract this 
background from the remaining fibers. Our algorithm did 
not make any a priori assumptions about the location of the 
forty sky fibers, but we confirmed after the fact that the se¬ 
lected fibers were indeed at the edges of the IFU, thus avoid¬ 
ing any galaxy light which could have otherwise affected our 
sky calculation. 
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Table 1. Comparison to Broad-Band Photometry 


Source U B V 


15.4±0.2 15.66±0.07 15.02i0.05 

15.5i0.1 15.6±0.1 15.1±0.1 


This Work 


Kim (2007 


Because the surface brightness in some regions of our 
target is more than 10 x fainter than the sky background, 
uncertainties in sky subtraction dominate the uncertainties 
in our flux measurements. Therefore, we adopt the standard 
deviation of the fluxes of the faintest forty fibers as the un¬ 
certainty in the sky background and propagate it forward in 
parallel with the flux data. 


2.5 Field-of-View Reconstruction 

To reconstruct the field-of-view at each wavelength we must 
take into account the relative position of each fiber with 
respect to the other fibers (provided by the McDonald Ob¬ 
servatory) as well as the sky position of the IFU in each 
exposure (provided by the guide camera astrometry). For 
the purposes of our project, we mapped the spectral infor¬ 
mation from every fiber in every dither into an x-y-X spectral 
cube with a 0”3x0”3 plate scale and a 3A wavelength scale. 
This plate scale oversamples the 2" resolution imposed by 
typical seeing at Mt. Locke. It is, nonetheless, quite appro¬ 
priate since a plate scale closer to 2" would badly under¬ 
sample some of the resolution elements because the fibers 
are circular (as opposed to square) and arranged in a hexag¬ 
onal configuration (as opposed to a square grid). Likewise, 
we have chosen a 3A wavelength scale, oversampling the 
5.1 A resolution by almost a factor of two. A choice closer 
to 5.1 A would effectively blur out by almost a factor of two 
narrow spectral features straddle pixel boundaries. 

We used a weight function for each fiber that represents 
the fraction of the exposure time during which that fiber il¬ 
luminated a given location within the image plane. For each 
exposure we determined the IFU position and drift radius 
from the guide camera images taken during that exposure. 
Our guide camera data as well as the VIRUS-P handbook 
indicate that the instrument drifts by roughly 0”5 during 
exposures. Each pixel in the reconstructed x-y-X spectral 
cube is a thus weighted average of all calibrated IFU spectra 
from all of the exposures of that target. The calibrated x- 
y-X spectral cubes are our primary data products, encoding 
the flux at every location within our targets at every wave¬ 
length covered by VIRUS-P. As described in Section |2.4[ 
the uncertainties in the x-y-X spectral cube are dominated 
by uncertainties in the sky subtraction, typically around 
a few xl 0 ~ 19 erg s - 1 cm _ 2 A _ 1 arcsec -2 . 


2.6 Spectrophotometric Calibration 

To correctly calibrate our target data we observed several 
standard stars at different airmasses each night. In order 
to ensure that the standard star flux is not truncated by 
the 4"16 fibers or the Vs AH factor of the VIRUS-P IFU, 
we compared the reference spectra to aperture photome¬ 
try from a reconstructed x-y-X spectral map of the field-of- 



Figure 2. Spectrum of region A' (see Section |3.2| l with bright 
emission lines marked; the other HII regions in UGC 628 show 
similar emission line characteristics, albeit at lower S/N. 


view. In cases where the standard star’s field-of-view was not 
completely covered, we reconstructed it using the azimuthal 
symmetry of the PSF. 

As a final step in our calibration, we corrected the fluxes 
in each wavelength bin in our x-y-X spectral maps for fore¬ 
ground galactic extinction using the Cardelli extinction law 
( |Cardelli et al.|1989|) and E( B-V)=0.0376, derived from the 
dust maps in Schlegel et al. (19981. 

As a check on our calibration, we passed the VIRUS-P 
data for UGC 628 through UBV filter transmission curves to 
recreate a UBV image. We compare the whole-galaxy UBV 
photometry to established values (Kim 2007) in Table |T] 
the photometry in all three bands matches to within 0.1 
magnitudes. 


3 ANALYSIS 


3.1 Emission Line Images 

In Figure [2] we present the spectrum from a region of 
UGC 628 particularly bright in line-emission. In order to 
identify regions with bright line emission and take advan¬ 
tage of the information available from line ratio diagnostics, 
we created emission-line-only slices from the x-y-X spectral 
cube centered on the [OlI]A3727, H 7 , H/?, [0III]A4959, and 
[O III] A5007 lines. Each slice is the sum of all the flux within 
an llA filter centered on each emission line minus half the 
average of the sums of two 22 A-wide continuum slices, each 
located 6.5A away from the emission line. The llA-wide 
filter generously includes all of the 5.1 A PSF, and the flank¬ 
ing continuum filters effectively take into account any local 
slope in the stellar continuum. Additionally, we account for 
the stellar absorption underneath the Balmer emission lines 
by adding 2 A of continuum flux back into the emission lines. 

This method of accounting for the stellar absorption 
underneath the emission lines was used by |Kuzio de Naray | 
et al. (2004) to correct H/l fluxes for the stellar absorption 


in a sample of six LSB galaxies and is motivated by the 
careful analysis of H II regions in HSB galaxies by |McCall| 
et al. (1985 I and Oey & Kennicutt (19931. Additionally, this 
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choice for our data is independently supported by the agree¬ 
ment between the H<5/H/3 and He/H/3 ratios (see Section 
3A\. 

To visually represent the distribution of emission¬ 
line luminosity with respect to the continuum emission 
of UGC 628, we present in Figure [ljffi) an H/3 image of 
UGC 628 with the stellar continuum overlayed as isophotes. 


3.2 Region Identification and Diffuse Emission 

Visual inspection of Figure [IJb) strikingly reveals HII re¬ 
gions arranged along the spiral arms of UGC 628. To analyze 
the sites of star formation as distinct entities we identified 
16 discrete regions by visual inspection, labeled A-K in Fig¬ 
ure |3]( a). These regions all met the criteria that they have 
centers 5er above the background, are roughly symmetric, 
and are separated from the next nearest region by a local 
minimum. 

Linking these regions to canonical HII regions is diffi¬ 
cult because the 2" seeing at Mt. Locke corresponds to a 
physical scale of 770pc at the distance of UGC 628. Com¬ 
plicating matters further, the physical size of HII regions is 
a thorny issue since no universally recognized method for 
delineating the edges of HII regions exists. Clearly, our data 
cannot resolve small HII regions such as the Orion Nebula, 
and the H/3 bright regions in Figure [3j a) should be thought 
of as HII region complexes and not individual HII regions. 
However, three regions, A, C, and L, have very bright cen¬ 
ters with sizes near our 770pc limit, which are much more 
likely to be genuine giant HII regions; we have identified 
these bright sub-regions as A', C ', and L'. The properties of 
these regions, described below, are summarized in Table [2j 
along with the diffuse H/3 emission and the net galaxy-wide 
H/3 emission. 


3.3 Extinction Maps 


To estimate the effects of extinction due to dust on the op¬ 
tical emission lines, we employed the H/3 and H 7 images 


and the extinction law described in Cardelli et al. (19891. 


We have restricted our extinction calculations to regions 
where both the H/3 and H 7 lines are 3<r above the threshold 


(see Section 2.51; this limits our analysis of the extinction in 


UGC 628 to only the emission line regions. 

We assumed an Rv of 3.1, as |CardeIi i et al.] (1 989j ) find 
this to be typical of the Milky Way, and an intrinsic H 7 /H /3 
ratio of 0.466, characteristic of Case B recombination with 
T = 10,000K ( |Brocklehurst|[l971| ). This information, along 
with our measured H 7 /H /3 map, allowed us to create Aa 
maps for each of our observed emission lines; as an example, 
we present the Ay map in Figure [3|b). 


3.4 Confirmation 


As described in Section |2.6| we confirmed the spectropho- 
tometric calibration by comparing our data to established 
broad-band photometry (Kim]j2007l; furthermore, with the 
emission-line images in hand, we performed two additional 
tests which confirm our calibration as well as our continuum 
subtraction. 

First, we compared the strength of the [O III] A4959 and 



Figure 3. (a, top left) H/3 only, with HII regions identified 

and labeled. Primed regions are particularly bright knots within 
larger regions, (b, top right) An Ay extinction map of UGC 628; 
note that the map is incomplete, as line emission is only detected 
from some areas, (c, bottom left) A log(0/H) metallicity map 
of UGC 628 derived using the R 23 method adopting the lower- 
branch value (see Section |3.5| )■ (d, bottom right) The comple¬ 
mentary upper branch log(0/H) metallicity map of UGC 628. 


[OIII]A5007 lines. The transition probabilities give a theo¬ 


retical ratio of [O III] A5007/[O III] A4959 = 2.98 (Storey & 
Zeippen|2000) and, since these lines are close in wavelength, 
this ratio is not strongly affected by reddening. We measure 
a ratio 3.11±0.04 in the 2" aperture centered on Region 
A' in the narrow-band images centered on [OlIl]A4959 and 
[O III] A5007. 

Second, the emission lines in Region A' are strong 
enough for us to detect the H<5 and, marginally, the He 
line, giving another set of fixed line ratios. Using the ex¬ 
tinction curve calculated from the H/3 and H 7 fluxes, we 
de-extinguished the H/3 flux measured in Region A ', trans¬ 
formed it to the appropriate values for H<5 and He using the 
Case B T = 10,000K ratios of 0.256 and 0.158, and finally 
re-extinguished them using the same extinction curve. The 
values we calculated in this fashion are 1.0 ± 0.3 x 10~ lb 
and 5.8 ± 0.2 x 10 - 1 'ergs _ 1 cm -2 for H<5 and He, respec¬ 
tively; the values we measured are 1.2 ± 0.4 x 10 -16 and 

9.4 ±4.3 x 10” 1 'ergs - 1 cm -2 . We treat this excellent agree¬ 
ment with the H<5 line and the acceptable agreement with 
the marginally detected He line as a confirmation of our data 
processing and choice of extinction curve. 


3.5 Optically Derived Properties 

Table [2] lists the properties of the emission-line regions. The 
H/3 luminosity was determined via extended source photom- 
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Table 2. Properties of HII Regions 


Region 

Rgai 1 

B-V 2 

log(L H/ 3) 3 

SFRh/5 

Ay 

l°g(0/H) 


(kpc) 

(mag) 

(erg s _1 ) 

(Mq yr- 1 ) 

(mag) 

lower 

upper 

A 

13.5 

0.641 

39.4 

0.10 ± 0.06 

0.902 

-4.0± 0.11 

-3.4 ± 0.040 

Al 

13.7 

0.740 

38.6 

0.026 ± 0.016 

1.19 

-4.0± 0.10 

-3.4 ± 0.038 

B 

11.4 

0.535 

38.7 

0.048 ± 0.030 

1.62 

-4.1± 0.15 

-3.3 ± 0.061 

C 

12.9 

0.587 

39.1 

0.062 ± 0.039 

1.03 

-4.1± 0.13 

-3.3 ± 0.052 

Cl 

10.7 

0.695 

38.6 

0.013 ± 0.0082 

0.740 

-4.1± 0.12 

-3.3 ± 0.047 

D 

9.28 


38.5 

0.0033 ± 0.0020 


-4.1± 0.066 

-3.3 ± 0.026 

E 

3.84 

0.605 

38.8 

0.016 ± 0.0099 

0.296 

-4.3± 0.11 

-3.1 ± 0.057 

F 

9.14 

-0.0855 

38.8 

0.016 ± 0.0010 

0.439 

-4.2± 0.12 

-3.2 ± 0.055 

G 

10.0 

0.618 

38.7 

0.0071 ± 0.0045 


-4.2± 0.11 

-3.2 ± 0.050 

H 

7.60 


38.7 

0.012 ± 0.0075 

0.376 

-4.3± 0.11 

-3.2 ± 0.049 

I 

7.08 


38.6 

0.036 ± 0.022 

1.55 

-4.2± 0.16 

-3.2 ± 0.067 

J 

12.5 


38.00 

0.0011 ± 0.00068 


-4.1± 0.12 

-3.3 ± 0.055 

K 

2.09 


38.6 

0.026 ± 0.017 

1.25 

-4.3± 0.15 

-3.1 ± 0.076 

L 

7.07 

0.513 

38.9 

0.023 ± 0.014 

0.582 

-4.1± 0.13 

-3.2 ± 0.059 

L / 

7.08 

0.463 

38.3 

0.0053 ± 0.0033 

0.527 

-4.2± 0.12 

-3.2 ± 0.054 

M 

15.7 

0.0965 

38.2 

0.0023 ± 0.0015 


-4.0± 0.071 

-3.4 ± 0.027 

whole 


0.64 

40.3 

0.45 ± 0.028 


-4.2i0.020 

-3.2 ± 0.011 

diffuse 



40.0 

0.17 





1 deprojected distance 

2 regions without B-V showed no statistically significant optical enhancement above neighboring pixels 

3 uncertainty dominated by 0.07mag uncertainty in photometric calibration 


etry in the H/3 image (see Section 3.1 >; the star-formation 
rates were calculated from the extinction corrected H/3 lu¬ 
minosity via the SFR/Ha in Kennicutt et al. (2009) with an 
Ha/H/3 of 2.87, again assuming Case B recombination with 
T = 10,000 K ( |Brocklehurst|[l97l I. 

The B-V colors of the emission line regions in Tabic [2] 
are the colors of the local optical enhancement in the x- 
y -A spectral data cube (referenced against Bessel B and V 
filters). The local B and V background for each region was 
calculated as the average value of the pixels falling on the 
perimeter of each region. All 16 H /3 bright regions show at 
least some some optical enhancement, typically much bluer 
than the surrounding areas. 

Table [2] also lists log(0/H) metallicities calculated 
via the R 23 quantity derived from the [Oll]A3727, H/3, 
[O III] A4959, and [O III] A5007 lines using the analytic formu¬ 
lae in [Kuzio de Naray et al.| (200 4). The R 23 quantity is a 
bi-valued function of log(0/H) metallicity; for low metallic¬ 
ities, increasing metallicity increases the strength of oxygen 
lines; for high metallicities, increasing metallicity decreases 
the strength of the oxygen lines due to increased cooling via 
infrared lines. 

Without the [NII]A6584 line to enable discrimination 
between the high metallicity ‘upper’ branch and the low- 
metallicity ‘lower’ branch, we present both metallicity values 
in Table [5] 


4 RESULTS 

4.1 Distribution of Star Formation 

The galaxy-wide star-formation rate measured via the ex¬ 
tinction corrected H/3 flux is 0.45 ± 0.028 Mgyr -1 . As dis¬ 


cussed in Section |3.5| we calculated this value from the 
extinction corrected H/3 luminosity using the case B re¬ 
combination Ha/H/3 ratio of 2.87 and the star-formation 


rate/Ha luminosity relationship from Kennicutt & Evans 


(2012). This star-formation rate is slightly above the aver¬ 
age of the LSB galaxies in [Kim ([2007J ) but we ll below the 
Milky Way’s 1.9 Mgyr -1 ( Kennicutt fe Evans||2012 1. Since 
Ha luminosity is the de facto standard when comparing star- 
formation rates between galaxies, we will make extensive 
comparisons between the Ha luminosities we derived from 
our H/3 fluxes and the Ha luminosities available for other 
galaxies in the literature. 

The most striking aspects of Figure [l]]b), Figure [3]] a), 
and Table [5] are 1) the distribution of star formation in 
UGC 628 into discrete clumps and 2) the alignment of these 
clumps with the bar and spiral arms. While our 770pc reso¬ 
lution limits our ability to resolve the smallest HII regions, 
we can say that the H/3 bright regions that we have identi¬ 
fied account for 50% of the observed H/3 luminosity but only 
15% of the area out to the same deprojected galactocentric 
distance (adopting the 56° inclination from Kim 2007)). 
Star formation in UGC 628 appears to be extremely con¬ 
centrated, though noticeably absent in the galactic center. 


Indeed, Region A, by itself, accounts for 20% of the star 
formation in this galaxy, and Region A' accounts for 20% of 
the star formation in Region A. Region A' is particularly 
interesting because it is concentrated even at the 2" scale 
testable by the VIRUS-P IFU (see Figure g a)). 

To place Region A' in context, we compare in Fig¬ 
ure |4jc) the calculated extinction-corrected Ha luminosity 
and the upper limit of the size of Region A' to HII regions in 
nearby galaxies (M33, M31, LMC, SMC, M82, M101, M51) 
as well as in the Milky Way (M42, M8, Carina, W49) drawn 
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Figure 4. (a, top left) An enlargement of Figure [3ja) high- 

lighting Region A!. Region A! is defined by the black 2" circle, 
which has been centered to optimize the flux captured from this 
somewhat lopsided distribution. Note that the H/3 enhancement 
is compact at the scale of 2 ", the seeing limit of Mt. Locke. (b, 
bottom left) An enlargement of the SDSS image in Figure Qa) 
identifying Region A! . (c, right) A comparison of the diameter 
and unobscured Ha luminosities of HII regions in nearby galaxies; 
marked on the plot is the calculated unobscured Ho luminosity of 
Region A ', the seeing-imposed upper limit on the size of Region 
A', and, to guide the eye, an Lhc* R 3 line. Since the Ho lumi¬ 
nosity of Region A' is in line with the upper limit on the diameter 
of Region A !, it is highly likely to be a genuine giant HII region. 


from previous works (Kennicutt 1984 Viallefond & Goss 
1986 van der Hulst et al.||1988 ). To guide the eye, we have 


also plotted an Lho oc R 3 line. 

The fact that so many of the HII regions from nearby 
galaxies fall near the Lh q oc R 3 line partly reflects an inher¬ 
ent size/luminosity relationship, but also reflects the fact 
that divisions between HII regions are often somewhat ar¬ 
bitrary. Many HII regions which can be classified as sin¬ 
gle structures can just as easily be subdivided into smaller 
structures, often down to the level of individual Stromgren 
spheres. The luminosity and size upper limit intersection of 
Region A' falls very close to this relationship, implying that 
it is likely either a single HII region with an unobscured 
log(LHa) ~ 39.6 or an agglomeration of several HII regions 
immediately adjacent to each other. Since the divisions be¬ 
tween immediately adjacent HII regions are fairly arbitrary, 
these two scenarios are not meaningfully different in terms 
of size and intrinsic Ha luminosity. 

Region A', then, is likely a giant HII region with a star- 
formation rate almost twice that of the largest HII regions in 
the Milky Way; while such a region appears extraordinary in 
an LSB galaxy with a global star-formation rate more than 


three times lower than that of the Milky Way (Kennicutt 


& Evans 20121, it is quite in line with the late-type (Sm) 


morphology of UGC 628. 

The link between the morphology of a galaxy and its 
HII region luminosity function is well addressed by |Caldwell| 
et al. ( 199l| , who carefully reconstruct the luminosity func¬ 


tions of HII regions in several nearby early-type (Sa) spirals 
and compare their findings to the similarly constructed lumi¬ 
nosity functions of late-type spirals described in Kennicutt 


et al. (1989). One of the major conclusions drawn by Cald- 
well et al.|(|1991|) is that star formation in early-type spirals is 


much more evenly distributed than in late-type spirals, and 


that the slope of the high end of the HII region luminosity 
function is much steeper in early-type spirals. Normalizing 
for differences in global star-formation rate, early-type spi¬ 
rals are much less likely to have bright star forming regions 
(log(LHa) ~ 39) than late-type spirals. 

The strong clustering of star formation in UGC 628, 
highlighted by the presence of Region A\ implies a shal¬ 
low, top-heavy HII region luminosity function consistent 
with other late-type spirals. While the argument is less com¬ 
pelling that Regions C' and L' are single HII regions, the 
presence of such highly concentrated star formation supports 
the conclusion that UGC 628 has a very shallow HII region 
luminosity function. 

The HII region luminosity function would seem, then, 
to be much more closely linked with morphological type than 
it is with global star-formation rate since the global star- 
formation rate of UGC 628 is quite low for a late-type spiral. 
This is in line with the findings of Schombert et al. (20131, 
where they construct an HII region luminosity function us¬ 
ing 429 HII regions in 54 LSB galaxies (mostly irregulars), 
and find that it matches the HII region luminosity func¬ 
tion for HSB irregulars presente d in |Youngblood fc Hunter| 
(19991. Indeed, Schombert et al. (20131 compare net galaxy 
Ha luminosities to the Ha luminosities of the brightest HII 
regions within each galaxy, and UGC 628 and Region A' 
would fit squarely within their data. 

That said, UGC 628 and the LSB galaxies in|Schombert| 


et al. (2013) both share one significant difference when com¬ 


pared to many HSB spirals: a lack of detectable star forma¬ 
tion in the core. Schombert et al. (2013) demonstrate with 
their sample of 429 HII regions in LSB galaxies that the cen¬ 
ters of LSB galaxies are generally devoid of HII regions and 
suggest that this may be an indicator that star formation in 
LSB galaxies is driven by local fluctuations in gas densities 
rather than global patterns. Nonetheless, star formation in 
UGC 628 is clearly influenced by some global process, as 
indicated by the strong alignment of emission-line regions 
with the bar and arms. Indeed, arms and bars appear to be 
the sites of star formation in LSB galaxies just as they are 
in HSB galaxies since |Kim ([2007) finds that nearly half of 
the spiral LSB galaxies in his sample show alignment of HII 
regions with bars and/or spiral arms. 


4.2 Upper vs. Lower Branch Metallicities 

While absolute metallicities cannot be directly determined 
using the R 23 strong-line method because of the ambiguity 
between the upper and lower branches, context clues favor 
the lower branch for UGC 628. For example, an examina¬ 
tion of Table [2] shows that the most compact and vigor¬ 
ous star forming regions, A and A', have the highest lower 
branch metallicities and the lowest upper branch metallici¬ 
ties. Since star forming regions are the sites of enrichment is 
highly unlikely that the most rapidly star forming regions in 
UGC 628 would also have the lowest metallicity. A similar 
argument can be made for Regions B , C, and C, the next 
most vigorous star forming regions in UGC 628, all of which 
have particularly high lower branch metallicities and low up¬ 
per branch metallicities. The upper branch interpretation of 
R 23 is essentially impossible to reconcile with the relative 
star-formation rates of the specific emission line regions. 

Additional context clues can be found by comparing 
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Figure 6. (a, left) A comparison of log(0/H) metallicity and Ay extinction for individual HII regions in UGC 628 and the LSB 

galaxies reported in |McGaugh| |T994[ ). Both upper and lower branch metallicities are plotted for UGC 628. For the HII regions from 
|McGaugh| < |1994| ), each host galaxy is coded with a distinct symbol; HII regions with identical symbols are from the same galaxy. Typical 
errorbars for UGC 628 are shown in the lower left, (b, right) Galaxy-averaged log(0/H) metallicity versus Ay extinction for UGC 628, 
Ha: selected galaxies from the KISS sample in |Young et al.| < |2014[ ), and LSB galaxies from |McGaugh| < |1994| . Both upper and lower branch 
metallicities are shown for UGC 628 (see Section |4.2| >. Since we only measure extinction in regions of UGC 628 where H/3 and H 7 are 
detected, the Ay values plotted are upper and lower limits. The upper limit is a flux-weighted average excluding regions where H/3 or 
H 7 were undetected, and the lower limit is a flux-weighted average including all areas where H/3 was detected, assuming zero extinction 
in the absence of H 7 . Similarly, the Ay values reported in |McGaugh| ( fl994[ ) are for HII regions, and should be taken as upper-limits on 
the galaxy-averaged extinction. 
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Figure 5. Galaxy-averaged log(0/H) metallicity versus B-band 
absolute magnitude relationship in UGC 628 (upper and lower 
branches) and a sample of LSB galaxies taken from |Kuzio de| 
|Naray et al.| ( |2004| l compared with that found in two samples of 
predominantly HSB galaxies: volume-limited Ha selected galaxies 
from the KISS proje ct [Young et al.||2014] , and spirals galaxies 
from [Zaritsky et al.| (| 1994 , as well as a diverse sample of galax¬ 
ies from Dcnico lo et al.| ~ 2002| l. Both upper and lower branch 
metallicities are shown for UGC 628 (see Section |3.5] . The up¬ 
per branch metallicity falls well within the locus of spiral galaxies 
and volume-limited galaxies, while the lower branch value is more 
in-line with other LSB galaxies and similar to the metallicities of 
several dwarf galaxies from|Denicol 6 et~al. (|2002[). 


the metallicity of UGC 628 to other galaxies. In Figure [5] 
we compare the metallicity versus B-band absolute magni¬ 
tude relationship in UGC 628 (upper and lower branches) 
and a sample of LSB galaxies taken from |Kuzio de Naray] 


et al. (2004) with that found in two samples of predom¬ 


inantly HSB galaxies: volume-limited Ha selected galaxies 
from the KISS project (Young et al. 20141, and spirals galax¬ 
ies from Zaritsky et al. (|1994 1, as well as a diverse sample 
of galaxies from |Denicolo et al.| ( |2002| ) . The well-established 
metallicity/luminosity relationship is clearly apparent, with 
LSB galaxies falling on the lower envelop. The galaxy-wide 
average upper branch metallicity of UGC 628, plotted as 
a red triangle, is similar to the metallicities of HSB spirals, 
while the lower branch metallicity, plotted as a blue triangle, 
is more consistent with other LSB spirals. 

A comparison of metallicity to dust extinction also fa¬ 
vors the lower branch interpretation. The relationship be¬ 
tween dust and star formation is complex: dust is an indica¬ 
tor of the dense molecular clouds that host star formation, 
it plays a role in the cooling processes that determine the 
IMF, and it is comprised of the metals that result from the 
cumulative star-formation history. In Figure |6](a) we plot 
log(0/H) metallicity versus Av extinction for the HII re¬ 
gions in UGC 628 for which we have both measurements as 
well as HII regions in other LSB galaxies described in|Mc-| 


Gaugh (19941. Again, both upper and lower branch metallic¬ 


ities are plotted for UGC 628. The data show a wide scatter, 
but there is a conspicuous absence of points in the lower- 
right corner, consistent with the idea that metals are re¬ 
quired for the formation of dust grains. 


The lower branch metallicities of UGC 628 fall within 
the primary locus of points in Figure [6j a), while the upper 
branch metallicities fall at the very upper extreme of HII re¬ 
gions in LSB galaxies. Indeed, if the upper branch metallici¬ 
ties were correct, that would imply that some HII regions in 
UGC 628 would have solar or super-solar metallicities with 
nearly zero extinction. 
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In Figure | 6 |b) we plot galaxy-averaged metallicity and 
extinction for UGC 628, the Ha selected galaxies from 
Young et al.| (|2014|), and the LSB galaxies from McGaugh 


(1994 1 . The metallicities and extinctions for the Ha selected 
galaxies were determined via long-slit spectroscopy, while 
the metallicities and extinctions for the LSB galaxies are an 
average of individual HII regions within each galaxy. Again, 
we represent both the upper and lower branch metallicities 
for UGC 628. The galaxy-averaged extinction for UGC 628 
was calculated in two ways. First, a direct H/3 flux-weighted 
average of the HII regions, and, second, an H/3 flux-weighted 
average assuming an extinction of zero for areas that are 
detectable in H/3 but too faint in H 7 for a proper extinc¬ 
tion calculation. The first method represents an upper limit 
on the true galaxy-wide average extinction since the regions 
which are brightest in emission lines are systematically more 
heavily extinguished. Likewise, the second method repre¬ 
sents a solid lower limit on the galaxy-wide average extinc¬ 
tion since the areas which are emission-line faint likely have 
at least some level of extinction, albeit systematically much 
lower than the emission-line bright regions. The true galaxy- 
averaged extinction falls between these two limits. Thus 
UGC 628 is represented as four points in this figure, two 
for the upper branch and two for the lower branch metallic¬ 
ities. The galaxy-averaged extinctions for the LSB galaxies 
from McGaugh (19941 have a similar ambiguity, however we 
are only able to plot the upper limits on extinction as data 
is only available for the individual HII regions. 

A correlation between metallicity and extinction is ap¬ 
parent among the primarily HSB KISS galaxies; the LSB 
galaxies also follow this correlation, but fall on the lower 
envelope. The lower branch metallicity points for UGC 628 
decidedly bracket the LSB galaxies. While the true galaxy- 
averaged extinctions of some of the LSB galaxies likely fall 
to the left of their plotted upper limits, the majority of these 
extinction values are likely not far from their plotted upper 
limits since the half-magnitude shift that would be needed to 
bring all of the extinction values in-line with the Ha selected 
galaxies would also lead to physically implausible negative 
extinction values for several galaxies. 

Based on the context presented, it is our conclusion that 
the lower branch metallicities are much more favored than 
the upper branch metallicities. 


4.3 Metallicity Gradient 

In Figure [3jc,d) we show metallicity maps for both lower 
and upper branches, and in Figure [7] we plot the average 
metallicity within each of the H/3 bright regions listed in 
Table [2] against their galactocentric radii. We find a mono¬ 
tonic slope in both branches; the upper branch has a slope 
of — 0.020 ± 0.002 dex/kpc and the lower branch has a slope 
of +0.021 ± 0.003 dex/kpc. The difference in signs for these 
slopes is related to the behavior of R 23 with metallicity, as 
discussed in Section |3.5| Nearly all low redshift HSB spi¬ 
rals show a negative metallicity gradient, typically ~ -0.05 
dex/kpc, with higher metallicities near the centers ( |Zaritsky| 
et al.|1994 van Zee et al. 1998 Ho et al.|2015 1. If we adopt 


the upper branch values for UGC 628, then this matches 
typical spiral galaxies quite well. If, however, we adopt the 
lower branch values, which are much more supported by the 
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Figure 7. log(0/H) metallicity of individual HII regions within 
UGC 628 versus deprojected galactocentric distance. On the same 
plot, U-V color of annuli at constant galactocentric distances. 
Both upper and lower branch metallicities are plotted (see Section 
|4.2| ; the upper branch shows a negative gradient and the lower 
branch shows a positive gradient while the U-V color shows a 
strong gradient of blue to red. 


context described in Section [4.2| UGC 628 is fairly unusual 
for a spiral galaxy. 

The prevailing explanation for the ubiquity of negative 
metallicity gradients is the inside-out hypothesis, wherein 
stars form earlier and faster in the central regions of galaxies 
(e.g., Goetz fe Koeppen| l9921. This hypothesis is supported 
by observations of steeper metallicity gradients in spirals out 
to z~2 ( Jones et al.|2 010 Yua n et al.|2011 1 , suggesting that 
metallicity gradients were steepest near the peak of cosmic 
star formation. 

There are several significant exceptions to the ubiquity 
of negative metallicity gradients which also support this ex¬ 
planation. Kewley et al. (2010 1 find that close pairs of galax¬ 
ies have shallower metallicity gradients, and conclude that 
gas tidally driven from the outskirts to the centers of in¬ 
teracting galaxies likely flattens their metallicity gradients. 
Similarly, Cresci et al. (2010) identify three z~3 galaxies 
with inverted metallicity gradients, and suggest that, at this 
early epoch, prior to the peak of star formation but dur¬ 
ing the era of mass assembly, primordial gas accretion onto 
galactic centers dominates. 

LSB galaxies may constitute another significant excep¬ 
tion to the ubiquity of negative metallicity gradients, and, 
thereby, provide another key insight into galaxy evolution. 


de Blok & van der Hulst (1998 1 measure the metallicity 


profiles of three LSB galaxies using the R 23 method used 
here and the FIVEL method described in[De Robcrti s et al.l 
(19871, and are unable to find a metallicity gradient. The 
metallicity profile of UGC 628, however, is constrained much 
more tightly and out to more than twice the galactocentric 
scale radius because of the larger number of data points in 
our sample. 

The distribution of star formation in UGC 628 hints 
that inside-out formation may not be an appropriate model 
for LSB galaxies. We we noted in Section |4.1| star formation 
in UGC 628 is decidedly absent in the core, with the large 
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emission-line regions located in the outer spiral arms. A sim¬ 
ple interpretation of the positive metallicity gradient is that 
UGC 628 has experienced outside-in formation, and that the 
outer disk is more metal rich due to a longer and/or more 
vigorous history of star formation. This explanation is sup¬ 
ported by circumstantial evidence from other LSB galaxies; 


for example, Galaz et al. (2006) and Matthews & Gallagher 


(19971 find inverted (red-to-blue) color gradients in 20% and 


50% of LSB galaxies, respectively. Radial migration of ma¬ 
terial certainly plays a role in determining the metallicity 
gradient, but it is unlikely to be a large part of the explana¬ 
tion for LSB galaxies since many LSB galaxies are in nearly 
solid-body rotation(e.g., de Blok et al.|[2001a|b de Blok & 
Bosma 2002), and are expected to be fairly inefficient at 


transporting material radially. If anything, they are likely 
to retain relic metallicity gradients should such a gradient 
ever develop. 

A complexity arises when we consider the clear, mono¬ 
tonic negative U-V color gradient exhibited by UGC 628 
(over-plotted in Figure [7]). The negative color gradients seen 
in most HSB spirals are another line of evidence for inside- 
out formation; the negative color gradient in UGC 628 seems 
to signal inside-out formation, contrary to the metallicity 
gradient. One possible explanation for this apparent paradox 
is that star formation was more vigorous and more centrally 
concentrated in the past, but has since slowed and moved 
to the outer disk. This is in-line with comparisons between 


models and broad-band photometry (e.g., Zackrisson et al. 


2005) and equivalent widths (e.g., Vorobyov et al. 20091, 


which conclude that LSB galaxies, while blue, are many 
Gyr old. In this scenario the positive present-day metal¬ 
licity gradient establishes an upper limit on the inside-out 
star-formation at earlier times, determined by the specifics 
of fSM mixing. From a practical standpoint, it is more likely 
that a rigorous measurement of the star-formation history 
of UGC 628 and other LSB galaxies with flat or positive 
metallicity profiles could help constrain ISM mixing mod¬ 
els. 

A significant line of evidence supporting this hypoth¬ 
esis comes from the fact that UGC 628 contains a bar, a 


fairly rare structure in LSB spirals (Mihos et al. 1997). The 


scarcity of bars in LSB galaxies is likely due to the fact that 
solid-body rotators are particularly stable against bar for¬ 
mation; however because of this stability, should a bar ever 
form, it would likely persist for many dynamical times. The 
bar in UGC 628 may be just such a relic since UGC 628 
is almost in solid-body rotation out to ~3 kpc (de Blok & 


|Bosma|2002[ ), where the bar ends. If this is the case, the bar 
in UGC 628 signals that this galaxy was once more suscepti¬ 
ble to global patterns than it is now, and may have exhibited 
inside-out formation at some point in the past. 

If star formation in typical LSB galaxies favors a model 
other than inside-out formation, this would conveniently ex¬ 
plain the reversed color gradients and flat metallicity profiles 
seen in other LSB galaxies. UGC 628, in this context, would 
be unusual for an LSB galaxy, having experienced a desta¬ 
bilizing event many dynamical times ago which left it with a 
relic bar and negative color gradient, but a metallicity pro¬ 
file that has been inverted by subsequent non-inside-out star 
formation, fn the upcoming series of MUSCEL papers our 
group will present robust, spatially resolved star-formation 


histories needed to flesh out or refute such a scenario for 
UGC 628 as well as other LSB galaxies in our sample. 

4.4 Metallicity as a Tracer of Chemical Evolution 


The evidence presented in Section |4.2| indicates that the 
lower branch metallicity is a more likely interpretation of our 
data. Throughout most of the disk of UGC 628 the metal¬ 
licity is close to log(0/H) ~ —4.2 (76 solar), with the most 
metal rich HII regions peaking around log(0/H) ~ —4.0 
(V 4 solar). While this value is typical for dwarf galaxies, 
low metallicities such as this are rare in massive galaxies. 
For example, in Figure [5] there are many HSB galaxies with 
metallicities at or below -4, but the majority of them are at 
least two absolute magnitudes less luminous than UGC 628. 
We can compare UGC 628 to HSB analogs by calculating 
what its metallicity would be if it followed the well-known 
mass-metallicity relationship. For the sake of example, if we 
adopt for UGC 628 log(M*/MQ) = 10.65 from Kimjj 20071 
and the mass-metallicity coefficients from |Tremonti et al.| 
12004), we predict log(0/H)= —2.9, more than a dex greater 


than the -4.2 that we report. Clearly, UGC 628 deviates sig¬ 
nificantly from the mass-metallicity relationship that holds 
so well for HSB galaxies. As massive yet metal poor disks, 
LSB galaxies such as UGC 628 can best be placed in the 
context of galaxy evolution by considering the physical pro¬ 
cesses driving the metallicity, metallicity gradient, and the 
formation of dust. 

The low metallicities of dwarf galaxies are generally at¬ 
tributed to their low escape velocities, which make them 
more susceptible to losing their ISM through supernova feed¬ 


back (e.g., Hunter & Thronson 1995 Thuan et al. 19991. 


This has the dual effect of removing metal-enriched gas and 
also extinguishing star formation. However, unlike dwarfs, 
UGC 628 resides in a deep potential well that makes ISM 
blowout from supernova feedback an unlikely explanation. 
For example, |DekeI fc Woo| ( [2003] ) find that supernova feed¬ 
back is no longer effective at evacuating the ISM from ha¬ 
los more massive than 3 x 10 10 Mg; the rotation curve of 
UGC 628 indicates a dynamical mass of 6.5 x 10 10 Mg within 
13.8 kpc. This key difference makes massive LSB galaxies 
particularly interesting for models of galaxy evolution be¬ 
cause, with ISM removal ruled out, the low metallicity of 
UGC 628 and other LSB disks must be explained through 
the history of star-formation and gas accretion. 

With this in mind, it is tempting to try to explain 
the low metallicities entirely through minimal past star- 
formation and think of UGC 628 and similar LSB galax¬ 
ies as unevolved analogs to high redshift galaxies. Despite 
the great difficulty in measuring optical emission lines at 
z > 2 and the large intrinsic galaxy-to-galaxy variation, the 
trend of decreasing metallicity with increasing redshift is 
seen across a wide range of environments and galaxy types 


(e.g., Savaglio et al. 2005 Erb et al. |2006 Rafelski et al. 


2012; Bay liss et al.|2014||Sanders et al.|2015| ). The standard 

interpretation here is an evolution of chemical enrichment 
over time. 

The analogy between high redshift galaxies and LSB 
galaxies is, however, imperfect at best, with a growing body 
of evidence suggesting that the stellar populations of LSB 
galaxies are quite old. For example, |Zackrisson et al.] (2005) 
study the star-formation histories of nine LSB galaxies using 
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broad-band photometry and find that, while precise histo¬ 
ries cannot be determined from only broad-band photome¬ 
try, ages younger than 1.5 Gyr are categorically ruled out. 
Meanwhile, Vorobyov et al. (2009) compare Ha equivalent 


widths from a sample of LSB galaxies to a simulated LSB 
disk, and find that the LSB galaxies must be at least 5-6 Gyr 
old. Vorobyov et al. (2009) also predict that a truly young 
LSB disk would exhibit strong radial fluctuations in metal- 
licity due to local enrichment; quite to the contrary, the 
metallicity profile of UGC 628 is quite smooth and some¬ 
what shallow. Indeed, the variation of only 0.3 dex across 
lOkpc, more than three scale lengths, categorically rules out 
scenarios relying entirely on recent enrichment and implies 
that most of the enrichment occurred many dynamical times 
ago. Clearly, LSB galaxies have undergone significant evolu¬ 
tion, albeit along a different track than HSB galaxies. 

A more complete view of chemical enrichment (e.g., 
Lilly et ah] 20131 treats the current metallicity as an equilib¬ 


rium between infalling gas, star formation, and stellar feed¬ 
back. Recently Peng et al. (2015) find that the chemical 
difference in the stellar populations of massive galaxies is 
best explained by a model where the ISM of star-forming 
galaxies is continuously diluted by infalling gas, while the 
ISM of quenched galaxies functions as a closed box. Once a 
galaxy’s supply of infalling gas is removed, the metallicity of 
newly-formed stars rises rapidly until the gas is exhausted; 
this is in line with the consistently higher stellar metallicities 
that P eng et ah] p015| find in quenched versus star-forming 
galaxies in their sample of 26,000 objects. 

UGC 628 likely represents the corollary — a galaxy 
where star formation slowed (or remained slow) despite a 
continued gas supply. The ultimate cause of the slow star 
formation remains unclear, but, regardless, if much of the 
stellar population of UGC 628 was produced many dynam¬ 
ical times ago then infalling pristine gas would have ample 
time to dilute the ISM metallicity and push UGC 628 off 
the mass-metallicity relationship. As discussed in Section 
|4.3| the positive metallicity gradient along side the negative 
color gradient supports the thinking that UGC 628 is quite 
old, and enough time has passed for the metallicity gradi¬ 
ent to be inverted by more recent (though possibly slower) 
outside-in star formation. 

This scenario is in line with Figure [6] where we see 
that, for a given metallicity, LSB galaxies typically have 
higher Av extinctions than KISS galaxies, which are pre¬ 
dominately HSB galaxies; the ISM gas in LSB galaxies may 
have been diluted by infalling gas, but nothing has removed 
or destroyed the dust grains. Interestingly, Gerritscn fc de| 


Blok (1999) hypothesize that, with a low metal content, LSB 


galaxies have difficulty forming molecular clouds and be¬ 
come trapped in a low metallicity state. If this hypothesis 
is correct, newly formed galaxies in the early universe may 
have been effectively racing against infalling gas; those that 
failed to produce enough metals to keep the process going 
became trapped in the LSB state. 


and gas-phase metallicity. Our results can be summarized 
as follows: 

1) Star formation in UGC 628 is highly concentrated 
into a small number of large HII regions, with very little 
star formation in the core. This tentatively supports the 
burst-and-quench model, where LSB galaxies are kept blue 
by brief but concentrated bursts of star formation which are 
unable to build up a substantial stellar populations. 

2) The presence of the largest HII regions rules out 
a steep HII region luminosity function, making UGC 628 
more consistent with other late-type galaxies even though 
the global star-formation rate is more typical of early-type 
galaxies. This implies that the mechanism which suppresses 
star formation in LSB galaxies arrests the process of HII 
region formation early on, and that, once an HII region be¬ 
gins to form, its growth proceeds as it normally would in a 
late-type galaxy. 

3) The lower-branch R 23 metallicity is the most likely 
interpretation of the observed emission-line intensities; a 
typical value is log(0/H)~-4.2. This metallicity makes 
UGC 628 more metal poor than HSB spirals with similar 
masses but more heavily extinguished in Ay than HSB spi¬ 
rals with similar metallicities. 

4) The metallicity shows a positive gradient across the 
galactic disk, with more metal rich areas nearer to the edge 
of the disk. This is at odds with a more typical color gra¬ 
dient (reddest in the center, bluer with larger radii), which 
implies that stellar populations near the center of UGC 628 
have experienced greater cumulative star formation. It is 
likely that UGC 628 did experience inside-out formation 
long ago, and that the current outside-in mode of star for¬ 
mation has since inverted the metallicity gradient. This sce¬ 
nario demands that the stellar population near the center of 
UGC 628 is older than the characteristic time scale of the 
current star formation. 

5) The low metallicity and high Ay/log(0/H) suggest 
that infalling pristine gas has diluted the metallicity without 
destroying the dust grains. This implies that much of the 
metal content in UGC 628 was produced at least several gas 
infall times ago. 

These lines of evidence suggest that UGC 628 and simi¬ 
lar LSB spirals may be the outcome of galactic ‘overfeeding’; 
that is, disks where continued star formation is stalled out 
by a low metallicity ISM created by a continuous supply 
of pristine gas when their initial burst of star formation is 
delayed, suppressed, or simply insufficient. This scenario, 
as well as others, will be rigorously tested in the upcoming 
series of papers where we will fit candidate star-formation 
histories against our optical spectra combined with UV and 
IR photometry for each region within UGC 628 and the 
remaining LSB galaxies in the MUSCEL sample. 
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5 SUMMARY AND CONCLUSIONS 

We have utilized the emission line luminosities and ratios 
in UGC 628, the first target of our MUSCEL program, to 
determine the magnitude and distribution of star formation 
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